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ABSTRACT 



The cofTiput at I onal rnethod developed in this theeie 
permits the cal cl?1 at i on or the aerodynamic performance of a 
i^ino of arbitrary olanform^ Both basic and additional lift 



st eady , i nvi sci d = i ncompressi bl e f 1 ow. The method oses a 
grid system of control points over the wing s-e.mi— span. The 
circulation ov^er the wing is considered variable with 
discrete values at the specified grid points. Finite 
difference equations are utilized to determine these 
discrete values. Control point i ndet er mi nac i es are 
eva 1 ua t ed an a I y t i c a 1 1 y - riat r i x i n vers i on i s r equ i r ed f or 
solution bv the method oresented. Details of the matrix 
technique are developed in Ref, 4, A brief summary of the 
principal computat i onal »^elations is included. No 

nurr»erical results are yet available but are expected during 
the n^xt phase of this research. 
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Linear Coordinates < 1 versus 

An g Li 1 ar Coord i n a t es ( 0 - © ) 

Field Element at Singular Point 

Computa.t i onal Grid in (<^,©) Plane 

Gamma Matrix and Vector Grids 
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The over a I 1 ree-earch of 
has a three-fold purpoee, 
instruct the student about 
govern the aer odynami cs of 



B A C K G R Q U N D 

which this thesis faros a 
Firstly, it is intended 
the fundamental equations 
wings in i ncompressi bl e f 



i ow , 



Secondly, it teaches the student how to convert the 
f ur damen t al equations into a form suitable for numerical 
calculation on the digital computer, and to carry out the 
complex progr amrni nq required for this purpose. Thirdly, 



once the first two objectives have been accomp 1 i shed , the 
computer program itself provides a useful pedagogical and 
design tool for illustrating the effects of various wing 
design paramenters on the final aerodynamic performance of 
a wing, 

Lcdr, J. L, pB.rk^ made a creditable start toward the 
first objective in his thesis IRef, ID, However, in the 
limited time available, he was unable to make significant 
progress toward the second or third objectives. 

It was evident early in this present analysis that 
Parks' computat i onal technique had some flaws. The 
ex pectat i on was that these flaws would be readily 
discovered and soon corrected and that objectives two and 
three would be rather quickly and easily completed. That 
has not been the case: 



S 



the investigation actually developed 



Et w r" I ( 



obetaci es 



e n c o u. P Ti e i“” ET Q a p? g ^ e*. l_ h 



ieri SB 



earn 



?veni:uaiiy overcome 



jrooress was nei tiner quick nor 



easy- In every case^ it found that the seerniPi: 
obstacle was really attributable to eoiTie conceptual error: 
Once the error was found, the obstacle disappeared and th? 



i n vest i gator s gained a deeper insight. Some of the major 
obstacles that were encountered and eventual I v' overco-ne are 
discussed in the following paragraphs. 

Confusion about the proper type of computat i onal grid 
to use, whether staggered or unstaggered , and the grid size- 
necessary for accuracy was a major factor. An 8 by 3 
unstaggered grid over the semi -span is now used, for good 



=a-son . 



Another point of confusion was how to deal numerically 
with the singularity that occurs in the governing integral 
equation. One option is to avoid the problem by staggering 
the grid of field points with respect to the grid of 
control points. The other option, which was chosen for 
thi-=. I n vest i gat i on , resolves the i ndet er mi na.cy by rigorous 
analysis snd as a consequence allows the employment of a 
simple unstaggered mesh. 

A third point of confusion concerned the evaluation of 
the partial derivatives of the circulation function. 
Should analytical or finite difference methods be used to 
represent these der i vat i ves7* It was found that analytical 



n 



•7- 



di -f -f erent 1 ati on is incorrect and that finite aifferences 
?T?u=t be used = 

Stiii another confusing point concerned the validity of 
r ep r esen ting the c i r c u i a 1 1 on f un c t i on b y a F our i er ser i e s . 
Uihile this pj^ocedure is widely advocated in the technical 
literature, it was found that a much simpler and clearer 
formulation can be obtained otherwise^ 

Boun d ar y c on d i 1 1 on s also c a u sed some c on f u s i on , 

Special conditions apply at the leading and trailing edges 
and at the wing tip and midspan. These quite compleK 
boundary conditions have now been fully and rigorously 
analyzed and incorporated into the formulation of the 
probl em. 



s evident from the foreocino discus? 



esearch has amounted 



c n a 1 



lo a major eoucati on in Dasic 
aer Ddynami OS and in numerical methods. In these respects 
it ha s b een a rich! y r ewar ding e>? p er i en c e . 

However, in view of the foregoing obstacles and 



problems, the 


time schedule 


has of 


course been 


great! y 


del ayed f r cm 


that which was 


initially 


anti Cl pated . 


Thus , 


there are no 


final numerical 


resul t s 


at this par 


t 1 cul ar 


scaqe OT cne 


i n ves t i Q a 1 1 on . 


There is 


also no mor 


e time 



available to this i n vest i gator . Hence, this final aspect 
will have to be completed by some subsequent i n vest i gator . 
Never thel ess , what the present effort has produced, to aid 
any subsequent worker, is a verv soohi st i cat ed and refined 



10 



riL*. 1 1_ 1 iTi^ i_ pj ljg s • hi 1 5 rn^trioQ r? 3. 3 r»ow 
that it can be ccnfidenTiiy Bxo^izii^d tQ produce 



0VOI V0Q t:o line dci nn 
^ r ei 1 an 1 0 



a r* i-j <=?. i_ L- u r a Ti 0 * i p? c?. i i~ 0 ■=■ li i. t : 



1 i 



i i . i i^^ ? WUULJS^ I i UIM 



jhiB anaivBis i= a con 1 1 nuax: i on of a previous stugv 

conoucceCi c*v^ uonn L_ » r arKs « Li eu’cenanc CoiTirriancier a (_»ni iieo 
States Naval Reserves CRsta 1 j» 

The concept o-f a continuous vertex sheet of v'ariaDle 
strength over the wing, trailing to in-finity aTt or the 
wing is utilized to determine ohe value o-f the ci rcu.l an i on 
-function at a finite number of control points on the wing. 
The vortex sheet strength is restricted by the Kutta 
condition at the trailing edge of the wing and by the 
requi r ement of no flow through the wing at the specif leo 
control points. The present analysis is restricceo to 
steacy, inviscid, i ncompr essi b 1 e flow about thin wings with 
straight or swept leading and trailing edges. 

ng aerodynami c s involves two fundamental aspects. 
F i r =. 1 5 ^ d esi red basic lift o i st r i but i on may ne sp ec i f i ed 

over the wing and the cor responai ng wing camber must be 
determined- Second , a f 1 at plate wing may be specified 
with the resultant adai ti on a 1 lift distribution to be 
calculated. Once the additional lift distribution is 
i^nown, other pertinent aerodynamic parameters may be 
Get er mined; notably, the slope of the wing lift curve and 
the location of the aerodynamic center- Also, when the 
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> 3 3 1 ci } » O r*. O 2 TZ X LJ n i 2 2 i" t O 1 3 ^ r" 2 D LI L. 1 O fl 3 d r~ 3 D O T1 fl r% Pi *_J W F? ;? 71 Pj 3 

: r;r r e3nond i nq induced draq can readily be cal CLil area = 



? Pi i 5 ci n d 1 V* ^ A ■=? Q e V e i o p *= 



niernf 



:ai cui a-re 



epanwiee and cbordwi ee pre33u.re di str i bu.t i one over the wing 
which are neceesas^y in designing the wing etrLictLire. 

Althou.gh Parke' cai CLii at i one did not achieve a 
eat 1 3^ actor y result, he pointed OL?t areas o-f concern -for 
possible errors and indicated a need tor -fLirther 

i nvest 1 gat i on . 

where Parks L?sed a staggered oesh between control 
points and -Field points on the wing to avcio the 
i ndetermi nacy enccu.Pft ered whenever a control point and 
■Field point coincide, the present analysis resolves the 
1 psdet er iT'i nacy and works with control points and fieio 
points super 1 nip osed Lip- on each other. Also. the series 

solLition employed by Parks to represent the c i r cl-1 at i on 
di str 1 but-i on has been discarded. Instead^ this 

d i St r i DLft 1 on is now represented more simply and directly oy 
the val Lies ot the circu.laticn at the -Fieio 

themsel ves= 



>Gi nts 
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A fl W i IdtiL-.'i'nti I T 

The wing georpietry considered in this devel opment 
consists of a wing pianform symmetric about midspan with 
straight leading and trailing edges containing no 
di scont i nui t i es , except at mid span where aiicwance is made 
for sweepback- The wing has no control SL?rfaces= A 
typio^.l wing is pictured in Fig. 1= with the geometric 
coordinates labeled. The wing is completely specified by 
thiree factors, trie aspect ratio^ taper . ano tr^e leacinQ 
edge sweeo. 



coordinates on tne wino are normalize 






semi -span length of one unit. Sweepback at the wing tip is 
Lamda<A)- Eta (T^) is the normalized spanwise coordinate. 
Parameter Tau(^) defines wing taper. Xi (^ > is the 

normalized chordwise distance from the leading edge. 
Si gma (cT) is a sign parameter , -^1 for the right wing and — 1 
for the left wing. The mean chord(c) is a function of 
aspect ratio(Hp^, that is. c ~ 2/AP. inis anaiv* = i = 

involves two coordinate transf crmat i ons . The 

transf ormati ons proceed from the (x,y)— plane to the 

non— di mensi onal ( -plane, then to the (<p ^ 9 ) angular 

coordinate olane- 
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Figure 1. Wing Coordinates 
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B. NUN— D Irib.NS I ONAL COURUlNAlES 

The “first change ot variables involves a tr 
r o rn c n e « v Diane into c n e ' "1 ^ a Pi e « 



'ma.D 1 * 



foliowinq rei at i onsh i ps . 



non— SLID sc r i pcea van ao i es 






± j. Lj 



p Q 1 n c s ; 



wniie variaoies wi' 



SLiDscr 1 pr 



■er 



control points- 



y “ 


; °y = 


— i \ 

\ B A / 


>'p 


= Y\ ■; (•) 

' ri • 

— 


( 3 B 2 ) 


L 


= ACT'tV 


(3b 3) 


'"Lp 


= AXp 


(3b 4) 


c = 


(2/AR) C(1 ^T) - 2rTri_3 


( 3 B 5 ) 


a 

LI 


= (2/AR> C(1 +T) - 2'r>|_^j 


(3.6) 


X = 




(3. 7) 


.-4 


~ '“P^P 


(3. S) 


gxl 


= (Vo-dY\_ 


(3= 9) 


dc 


= - <2/AR) 2T'(Tdry^ 


(3= 10) 




= b X c d 1 + ^ d c 


(3. 1 i ) 




= A<rdv\ -i- cd\ -J- 1 r- (2 /hR) 2T'<rd Y\_ j 


(3. 12) 




L h — (2/hR) 2 'T' 1 J 


(3. 13) 


dx 


= cd 1 + 


(3b 14) 


r = 


(i — taper ratio)/ (1 + taper ratio) 


(3. 15) 


Thi 


s transTor mat i on simpii-fies the algebra 


or the 


obi em 


significantly. The chordwise dimension now 


var i es 



•from 0 at the leading edge to -*-l at the trailing edge, 
while the spanwise dimension varies -from —1 at the left 
wing tip to 1 at the right wing tip. 
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k ifD) 



\ o- - i / 



Lj 3 Hi'vLjULHH! L_-lJUr\iJ i r--lfn ! JZ.w* 

It i = desi^aole to trans-form ths wing pianform from -chs 
<|=>^)“pians into the (0=0)-piane» Tnis is accornp 1 1 shed 
Lising the rsi at i onshi ps 

^ = O - 5 -fr <1 — c os ) ; d I = 0,5 sin (t> d 0 

*'0^ = cos 0 ; = -sin 0 dd 

From equation (3,16). as \ varies -from zero ac the leading 
edge to 1 at the trailing edge, will '-^ary from zero to 
r espect 1 vel y - From equation <3-17), as r\^varies from “i at 

the left wing tip to at the right wing tip* 0 will vary' 

froim TT to zero, respect i vel y , These limits are important 
in subsequentiy developed integral equations. An ad dec- 
bonus of this t r ansf ormat i on is that it results in mesn 



points near the leadino and rr-aiiino eoaes ano 



VSr i ! ! U 



ICS to 



be closely spaced. 



owi nq 



n er resol u 1 1 on 



critical 


a r e ci s of a ns o. i y' r? i s - 


Fi g - 


2 shows 


the 


r e 1 at 1 on sh i p 


between 


the linsB.r coordir 


i^t es 




and 


the anqul ar 



coordinates ( 0,0) 



i — 




0^0 



0^'ir 



Figure 2. Linear Coordinates (^ ,^) versus 
Angular Coordinates ( 0 , 9 ) 
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Trailing ->*» Edge 



i V . bOvtnRNINH h.LjlJAfIQNb 

The vertical velocity iridLia^d at any particular control 
point consiete o-f that contribution made by the vorticity 
distribution over the wing piantcrm plus the contribution 
made by the trailing vortices behind the wing» Parks 
CRef. 1 j provides a quite detailed derivation ot the two 
general equations tor induced velocity. Both equations may 
also be round in many other aeronautical rererences such as 
Bertin and Smith [Ref. 2 j and Kuethe and Chow [Ref, 31, 
The fallowing form of these equations was utilized in this 



a n a i 'v ^ ^ ^ 



K 



k/j t k' ^“1 



r 1 [ 1 


(4 


■ dn 


> (n-if)! 


4 J 





wner e 






- v'p' 






and r is obtained by replacing x with in equation 

(4,3) . 



Equation <4.1) is basically the same as Park s equation 
(3.6) where ' is the velocitv induced at a control ooint 
due to the vorticity distribution over the wing surf ace » 
Equation (4.2) is the same as Park's equation (3-S) where 
w„-‘ IS the velocity induoBd at a control point due to the 
trailing vorticity behind the wing. Subscript t here 
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Winn, 



r 0 -f er 5 t. o V 0 1 u. 0 5 along t: n 0 i: r a. i i i n a 0 d g 0 o -f t: n 0 w 1 n g = It: 
10 d0sirabl0 to transtorm those eqLiations to equivalent 
f or n s in the a n □ u 1 a r c oar d i n a 1 0 s y s t em . 



A. CHANGE QF TO WING CQGRDINATES 

Circulation has the dimensions of veiociiiy times 
distance- Consider two separate circulation functions, one 
non-di mens! onal with respect to semi -soan and fi^ee srrearn 



veiocitv. ana one non— a 1 msnsi ona i wi tn resoeci 



10 (Tie ar- 



ch or d and free 

reT at i onsh i ps can i 

A 



stream velocitv- I he followinq. 



:nen be esrabiisned 






C/S = Z/Ht\ 



(4. 4) 



A 




Pil.n.) 






1 3 / 



t ^ 3 'tl. } 



iz-ubsiii nut 1 ng equations Ci. 



1 ) an d ( 3 3 i 4 ) into (4=6) q 1 ves 




( 



/ ) 



Kegrouping terms 




( 4 . 



The coefficients of d^ 



\nd of dr^must separately vanisr 



2u 



G 1 V' 1 n G G W •_» £ Q 'Li 3. G 1 Q n =■ W G 1 C Pi P G LI C S LG 



bxy 



A 1 /il! 

AR C Ul; 



z 

AR 






(4. V> 



I 4 = 1 i 



O) 



Now, the -f'Gil owl n-G e>?pression that occurs in eqLiation \4=1> 



'"P'Uxy 



( V ~ V 



.Ar 



'■p' 



( ^ = I 1 ; 



can oe exoressea in rerms of J a 



nd Y(^ as 



T Q1 i ows 






iS-fd 



C 4 = IZ) 



■Duping terfns in equiation (4, 12) and setting 
F = <>( - Kp) - <>l-7|_p)/^cr 
dx dy = c d^ drj^ 



<4. 13) 






hqjuation <4,i) may now be written in the following form 

f i . I 






±IA\ 

4trlARjy, 4 






lo; 



B. CHANGE QF P JO ANGULAR COORDINATES 

Utilizing eqLiations <3=16) and (3-17) together with a 
non— di mensi onai circu.iation fun-ction in the (^=0)— plane the 



■oi lowing tr an sf or mac I on can oe mace 






(4.16) 



( 



AL 

nl 









n = 




(4. 17) 




-sin 9 off) 





[quation (4.15) now becomes 



o tr 



w.. 



4.ttIar 



er-rr 






(4.23; 



Upon reversing the limits of integration and reversing the 
algeoraic sign of the integrand, equation (4.23) can be 



wr i t ten as 



rTTr'^ 



■p 




4 rrlARyy j r 

o 6 






Equat 1 on 


! A A \ 

\ — r a / A 


s the governing equ 


ac 1 on 


for che 


1 nduced 


vel oci ty 


at each 


control point due t 


G the 


summat 1 on 


of tihe 


effects o 


f every 


field point on the 


wi ng 


pi anf orm a 


Not e 


that al 1 


f a.c t or s 


inside the integral 


are a 


f uncT 1 on 


of winq 


p 1 anf orm 


geometry 


except for the two 


parti 


al den vat 


1 ves of 


the circu 


1 at i on f 


unction. Also note 


that 


equat i on 


( 4 . 24 ) 


becomes 


i n determ 


in ate when x = Xj_, 


an G 


> 

II 

r 


Thi s 



condition occLirs when the controi point and the -?ield point 
coincide and ie diecLieeed later^ 



C, TRAILING VGRTICITY TRANSFORNAT IGN 

Many o-f equations (3el) through (3.14) are identical at 
the trailing edge and simply require the subscript "t' to 
be added to the appropriate variable- However ^ since J is 
constant at -^1 along the trailing edge a tew ot the 
equations change slightly and are rewritten here tor 



Cl an ry 



(3-7) becomes = x^ + c 



(3.3) becomes Xp = Xj_^ Cp 



(4 - > 

(4-26) 

(3-12) becomes dx^ = L A - (2/ AR) 2TD (Tdr\^ (4-27) 

(3.14) becomes dx^ = ^CTdi^ (4-28) 

The circulation along the trailing edge ot the wing is 
a function of the spanwise coordinate only. Consider 



] ne 



de'/^i ap men t: is similar -co tnar pressnrsd aDove r or tns 

circulation the VJinq and results in the following 



•qo 



(T <y) = and c/s = 2/AR 

rv ^ y ^ ^ 2 /hr ) (v\) 



(4« 2V) 
(4=30) 



f 4 i)a. 


2 


/dR) 


Uh/ ' 


AR 


(da/ 


fii] - 


2 


/<)R\ 


Uh / 


AR 


Ur^J 



gV^ and dy = dr^ 



i ^ = o 1 ; 



^ 4 » -j.2 ; 



Tt ^'"'t "'P ^ 



< 4 . 33 



G = 1 - C (>i^ - > / r . 1 

r p r 



(4.34) 



Equation (4.2) may now be written as 

.+1 



I f2 






G /df^ 



4'fr lAR/y^ (r|_- VdY\> 



(4=35) 



Now, tr ansf ormi ng to angular coordinates 



(Y^) = Fl (6) 



(1^)°'^ ' (le] 



H . r.p; 



i 4 . O / ) 



Equation (4.35) now becomes 

rO 



-p" = 



4ir Iar) 4 K'^pl ^d© 



d 4 l 



(4.38) 
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Again, reversing the limits or integration ana the 
algebraic sign of the integrand, the integral equation ror 
the induced velocity due to the trailing vorticity becomes 

= — - — /— J / (4.39) 





\ ^ 1 




\k?sl 


Vo * 


[69j 



D. TOTAL VORTICITY EFFECT 

The total indL?ced velocity at a particular control 
point is the summation or the effect produced by the 

vorticity over the wing and the effect produced by the 

trailing vorticity- This means that for each control 
point 5 the right hand side of equation (4.24) must be 

evaiuatea at ever 
and summed, Thi= 
vorticity over the wing. 



field point over 


the entire wing surr 


ace 


um will give the 


effect p r od uc ed n y 


r h 


wing. Also, for 


each control poini:, 


the 


equation (4,39) 


must bs f^valuatsd 


at 


the entire trai 1 i 


ng edge of the wing 


and 



summed. This sum will give the effect produced by tne 
trailing varticity= Uihen both of these are added. tne 
total vorticity effect at that particular control point is 
evaluated. This may be expressed as a combination of 
equations (4,24) and (4,39) as follows, where Wp is the 
total induced velocity 



w = w 

■p -p 






( 4 . 40 ) 






■ iJi I 



INUti \ hiKiilNACi 



Equar i 


ons (4, 


24) and 


(4,39) may be evaluated 


by seven's 1 


Qi f f erent 


methods 


■ } n e 


relationship between the 


g earns tr i c 


pi ac ement 


of the 


con t r o 1 


points and -field points 


is pernaps 



Qne or the iTfOst important tactore- wing evmmetry about 
mi dspan permite the anaiyeie ro utilize control points ovor 
only one semi— span or the wing. The etfect of field points 
an d t r a i 1 i n g •v'or 1 1 c i t y .must i nc 1 ud e the entire wi ng 
pianform: however, symmetry B.baL\t micspan again permits 
utilization o+ numerical values over only one semi -span. 
The right semi -span is chosen for this analysis. 

As mentioned befoj^e:, equation (4.24) becomes 
i ndeter mi nate if the control point mesh and cne field point 
Tissh coincide. Equation (4,39) also becomes i ndet ermi nate 
when the span wise coordinate of the trailing vorticity and 
the concej^ned control point coincide. Note cnat G in 
equation (4.34) approaches zero as the spanwise cooraina.te 
of the control point and the trailing vorticitv converge^ 
One method to avoid these i ndetermi nat e points is to 
stagger the different meshes so that they never coinciae. 
Another method is to determine an analytical lim.it or 
principal value of the integral as the i .ndet .mi nat e poinc 
IS approached. The second method is chosen for this 
analysis and is p.f^esented in the following paragraphs. 
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INDETERMINATE POINT OVER NIN^ 



a 



? 11 r o 1 p o j 



tieiG DDinr 






numerator and the Genomi nator of eauation ( 4 =:^ 4 ) are bo~h 
z ero reeul x. i nq i n an i ndeter mi nacy . Si nee an i ndet ermi nacy 



a n i V o c c u r a -c 



c on e r o 1 o o i n r . 



confined to the right aarni“span= A] 
and sigma, the sign parameter » have 



Ml occurrences are 
g e o m e o r i c v' a r i a o i e s 
P o s i t i V e V a i Li e o n 



this section of the wine. bio, 3 deoicts an expioc 



of a mesn poir 
circLilation aT 



>inaularitVa I ht 



1 1 serene v a i Lie 



:hs singu.l ar point mList be determinea = 



conn nuoLis 



Recall that the circu.lation function itself i 
over the wing; therefore, it is continuous in Tne proximitv' 

V' a r 1 a D i e S- li a r? d v ^ }•” ^ 



of cne sincLiiar oox 



incf'oducea to represent small deviations i 
vicinity of the singularity. If AR 
contribution of the singLilar element to 
velocity at the singular points the fcliowin 



s in ^ a n G ^ = i r = c r? in* 



r ep resen cs 



me indLiceo 
r el at i onsh i p 



is established 







w n e r t* s y m o o 1 


^ i 5 an 


abbrevi at i on 


"c n a c 


stanos 


integrand of 


eoLiat i on 


( 4 , 24 > . Now , 


1 et 




Li = ( - 











for the 






Z/ 
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Figure 3. Field Element at Singular Point 



V = ( 9 - d^,> 

^ i ij. j V ) = F ■ / r ' 

where F' represents i 
(4=24) =. 



rne e>ipression i 



! > (_• ] 



3 ) 



ackens in eGustion 



H.n analysis or equ.a.TiiDn in ‘cne i iT.niedi ar e 

vicinity of the singularity results in the following 

K = (2/AR) L (i -j- T ) — 2 ti^J 1 (5=5) 



Kp. = Aq,-, + (2 /Ak)C( 1 T ) - 2^qpjlp 



\ Zi . ; 



<K-Mp) = A0\-qp>-(2/AR)C<l-^T) (i-lp)-2T(T^]-np?p) j (5.7) 



t\ = COS 6 = cos(8 ^v) = cos & - sin 0, 



V n a .} 



Id = 



(5-9) 



- no' = ■ ®£ 



(5- 10 



\ — 0.5(1 - CDS 0 ) = 0 . 5 C 1 — 



cos ) J 



= 0-5(1 - cos 0p -r sin 0p l? ) 



i L } 



)d = o.5(i - CDS 0 : 



(5= 12) 



( ] - Ip) = 0. 5 sin (j6p u 



VO- iZ} 



~ <lp+0. 5sin(^^ u)-y\plp (3,14) 



^p^p 0.551 n^p u - 5in0p v 

- sin 6p V 0.5 sin 0p u - r^plp 



<5.15) 



Simplifying and neglecting higher order terms in = 



•9 



and uv 15 



^'n}-Tlplp'^ = L sin (^„JU + (-sin 0p ) v + = ..(5.-16) 

SLibat i tuti ng equ.atians (5, 10) ^ (5c 13) , and (5= 16) into 
( 5 » 7 ) , 



( ~ M ^ ) — Alj -j- B V 



t. D = i / } 



^hers A and B are known functions of A? 6 n^ ^ ^P' " 
and AR. The denominator of the integrand of equation 
(4c 24) = in the immediate vicinity of the singular point may 
now be e?^ pressed as 



r‘- = L <Au -j- 



^ Bv)"^ + (—sin 



( 5 c IH) 



dnce is fixedc let —sin = C, then, 



T ^ -J I=- 

= r (Au + Bv)"^ + 



(5c 19) 



Assume, by analogy with the above development, that the 
numerator of equation (4c 24) may also be expressed as a 
function of u and v, namely 

F ^ = D u “^ £ V ( 5 c 2 0 ) 

T'lcw , switching to polar coordinates with ' the polar 
radius and ' c>i ' the angle, let u = ^ cos 04 arid v = ^ sin oC > 

The effect of the singular field element must be 
determinedc An integration over the entire element 



produces the required result 



The 



1 nteqr at 1 on 



covers 



total angle of 2TT and a polar rs.dius- varying from zero at 
the center of the element to its outer bounds * Equatione- 
•' 5 « 1 9 ) and iS. 20 ) ex pres sed i n p o 1 ar coord i n a t es become 

r " = L (Acqs oc — Bsi n oc ) si j ^ ~ (5=21^ 

= J (<^> p ^ (5= 22) 

where J (<^) is an abbreviated form of the term in brackets* 



h ' = L D cos oc + E sino63 



= H(^) ^ 






where H is also an abbreviated form* 

Note the palBir symmetry of each of the above functions, 
that is 

. J ( oc -f- It ) = J (oL) (c*. 25) 



H ( + 'Tt ) = — H 



Therefore, if AR' is the effect of the 




The value of this integral is zero 
pal^r symmetry. This proves that 
singular field element on the control p 
be neglected- 



i n deter mi nacy 







(O. 


27) 


by 


reason 


of 


odd 


the 


effect 


of 


the 


oi nt 


is zero 


and 


may 
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&- INE)ETERMIN.ATE POINT ALONG TRAILING EDGE 

Ab appr Qach&5 in equation (4»39) , G apprcaonee 

zero; therefore^ the integrand is i ndetermi nate . Note that 
(><• — K' ) does not aoproach zero. since is the 

T p ' • 

>?-ccordi nate or the center ct the concernec control point 
and is the k— coordinate at the trailing edge of the 

wing. In the i iHiTiedi at e vicinity of the singularity^ 






ecoiTies verv verv sfT^ai i , but remains finite; 



therefore^ it can be represented as a small value times a 
constant. Let this constant be (k- - therefore. as 

approaches zero, the following rel at i onshi ps may 
be established 



= e^(x^ - 5 where 



'H " 'V'"" ^ '’H. - np>" 
) - ( 1 ^ 



i 5 . LV) 



r^ = ^ (1 ^ 

= ( X ^ — X ^ ) ( 1 + e"^ / 2 . = . H . O „ T - ) 



( u . o U ; 



The integrand of equation (4.39) may now be expressed as 



(V'IpI 



X-t-X 



? - 



Xi-Xf 



R i ^ (xr L (X-t“Xp)(l i- 



(5.31) 



which can be reduced to 



(Xt \ I + 



(5.32) 



Thus, in the limit, as apprcachee ^ approaches 
zeror and the value of equation <5-32) becomes zero* This 
proves that the trailing vorticity contribution dir'^c.riy 
downstream ot a control point makes no contribution to the 
velocity induced at the control point. 



!-’KtLIMINAh;Y PRuyLEM ANALYsihi 



^onsidsr a oivsn winq DlanYortn tor wnich the three 



w. • cl ! • I 



known vai uas, Tne cor r sspona i nq 



M i i T •- 



= -r T }^ • J 






quani: lilies ra as aecer -tii nsc , hi 
mi: eg rands of equations (4,24) and (4,39) ar^ 
-for the partial derivatives ot the circulatic 
A mesh size adequate to repr"esent the 
bution. while not requiring an excess: 



u.anci Cl es 



encire wino soan is arhicrariiv' selei_ceG aPiO should aroouce 



me cf 



SViTnT-eCr'' 



accur acy , 



iS previously mencionec 



lep 


are 


on 


and 


ng 


are 


n 


the 


BK 


cent 


1 on 


» 


id 


X J. ! 


?unt 


Of 


1 n g 


nf 


‘SS 


the 


pro 


dues 




wi no 


, ^ 


made 



w 1 c n t: n e n umer i c a j 



uues on the rioht sefni*~span. 



A= wirJG GRID PQINTS 

This analysis involves /napping points on the wing rroo 
the physical (x,y) —plane to cor r espondi no points in the 
(0,6)— plane and vice versa. Fig, 4 depicts a uniform 
rectangular grid o-f points in the (^,0)— plane, constructed 
as ■j'Oiiows 

0. = I A0 = 26 = TT/S; i = 1,2,..., 8 (6.1) 



Y.4 




Figure 4. Computational Grid in (^,0) Plane 
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Trailing Edge 



0, 

_l 




■■ Ae 


< j-0, 


, 5) 




A6 


= 2cT = TT. 


'id: J = 1,2,=,., S 


(6 . 


'y ) 


r—i ^ •__! — 


al 


on g 


the 


1 e 


ad 


1 ng 


5- n Q "c 5 1 


.ling edges of the 


Wl ng 






( j 


) = 






( j 


) = 


TT 




\ Cj = 


3 ) 


wh i 1 e 


a 1 


eng 


the 


Wl 


ng 


1 1 


p and at 


mi dspan 










\ — 

* — 




e. 


(i 


) — 


Tr/2 




\ d* ■ 


4) 



jH 



The above equations will produce 64 grid point 
right semi -span with each point centered in its respective 



• lie 



j' subscript vanes span wise 

i ' subscript 

varies chord wise from leading edge to trailing edge. The 
unknown values of the circulation function at each of these 



rectangular mesh b i oc rc . 
from right wing tip to mid span while the 



For every field element on the right wing there exists 



a cor respond no mirror i me 



finer e exists 
ield element on the left 



m 



C\ 3 ± .- x 



ano 



itities witn 



wing. By utilising the symmeti^y between each field element 
and its mirror image, the range of the 'j' subscript can be 
restricted to B vice the 16 it would otherwise be. Fig, 5 
depicts the planform distribution of the circulation 
function with appropriate subscripts for 
cor r espondi ng vector represent at i on = All qi 
numerical values at respective grid points may be 
designated similarly. Should it become necessary to 
transform from matrix notation to vector notation or 
vice versa, the following r el at i ons»h i ps are used, wit.h 'k' 
denoting the .r-espective vector subscript- 

k = S(i — 1) j where ( k ( J ^ i ) (6 = 5) 



Lea^iitig Edge Leading Edge 



Righ-t Wiagtip 



r.(U) 


L(l>2) 


T(1,3) 


* « • 


rji,8) 












L(3.l) 


• 

* 

* 


rjin) 






» 










r.(84) 








r.(8,8) 



Midspan. 



Right Wingtip 



rv(i) 


L(9) 


r,(i7) 




R(57) 


au) 


L (10) 








L(3) 


s 

4 

4 








« 

* 










L(8) 








R,(6'7)- 



Mid span 



Figure 5* Gamma Matrix and Vector Grids 
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Trailing Edge ■ Trailing Edge 



6 ) 



I = k / H 1 < k / H by in t eq division) 









i) = f[.(k) 



(6.7) 



S. VARIABLE MGTATIGN 

The SL-becript notation introduced in the aoove 



paragrapns is used expensively in cde suGsequent: 



0. 1 ! ^ A V “ 



The foi lowing notation 
derivatives is al so used. 



f or var i ous f un c t i on s 

Note the different exoress3 



that czarr^Bpond to the same numerical value. 



The first derivatives of the cij^culation function 





( 





1 0 ) 



( 6 , 



i 1 ) 



The second derivatives at the leading and trailing edge 




( 6 . 12 ) 



(6. 13) 



3S 



The lift di str i but 1 on is designated as toilow= 
AC^) . , = P( 1 .1 ) = p..(k) 

p J , i * ^ V 



AC_ j “ Pi ^ J ' 

p _» L- 



( 6 . 14 ) 



saas 



tj 



t 



( 6 . 15 ) 

( 6 . 16 ) 



C= LIFT DISTRIBUTION 

Consider BernoL-lli's eqsuation tor t.he pressLs.re above 
the wing and below the wing 



Pn ^ <C 



/ ^ j r i j _i_ 
/ ^ y s_ V 



OO 



a , . / ~i ^ 



= D 



L 



+ <^/2)Cv^ - '!T>j/2j2 



i 6 . 1 / ) 



Rearranging and simp! i tying gives 



tevi. 



* i 



‘ fer, 



( 6 , I i 



A 



_ . -L 1 /iC 
J AR c y 



i. 6 . i ' 



and , 



^ IK] 

5 *irv 0 / 



^ O - xlV 



Now. combining terms from above 



AC 



lAR/c \6ir>. AR cl5in0/U!> 
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Ur, in the notation inx:rodLicea in section b 



P ( j 



(S\ 




) 1 


IarJ 


i CA) -^irv 0i J 



(b= 2:^) 



LEADING AND TRAILING EDGE REQUIREMENTS 

The ba\j.nds.ry conditions pertaining to equation (6-22) 
tor the basic or ideal ii + t, and the additional lift, 
determine the nature of the circulation function and its 
derivatives. The Kutta condition fixes the oressure 



coefficient at zero at the trailing edge. For 



one 



cionai iii-t. -rne oressure coe-rriciem 



inn nice at 



the leading edge. However, for the basic lift, 






pressure coefficient at the leading edge is finite, 
specifications together with equation (6.21) or (6=22> 
determine the nature of the circulation function. 

First, consider the Kutta condition which requires that 
the pv- assure coefficient be zero at the trailing edge. The 
Gsnomi nat or of equation (6.22) is also zero at che trailing 
edge, since ^ = TT and consequen t 1 y , sin ^ = O. This means 
tnat must not be finite or the quotient becomes 

infinite. Therefore, at the trailing edge must be zero, 

resulting in an i ndeter mi nacy . Using L'Hopitai 's rule to 
determine the limit, which must be zero, the denominatior 
becomes cos TT = — 1 and the numerator must be zero. This 
fixes the value of the second derivative of the circulation 
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t un c t i on with r e=p ec t to ^ a t 



?ro aiSQ, }Q Bumi'narize xi 



BQLiB.iiiiDn rorni 



Ac^} ^ < j) =0 















(6= v:^) 



i t. = ; 



’. a ^Z 



Next. consider the condition that the pressL-re 

coet-ficient is infinite at the leading eoge -for the 
additional lift. At the leading edge, 0 = O and sin = On 
therefore, the denofninator of equation (6.22) 

, is fini te a 



± =• ci •_! -rS i. ! i 



zero. Thu 
addi t i onal 1 i -f t . 

ACj3>i_ = P^Cj) 



the leading edge for 



^r-u.mmar i : 



= -f- OO 



(M ‘ ■ '' 



k <3 . 



( 6 . 2 /) 



Finally, -for the basic ii-ft the pressure coefficient is 
finite at the leading edge. A similar analysis to zhat for 
the trailing edge using L'Hopital s rule results in tre 
second dBri'/BtivB of the circulation function with respect 
to being a finite value. Summarized in equation form 



A( 



Cp)L = Pl(j) = (a/AR) 






( b - z: V / 



KtQU I RhIiEN I S 

tipj denoted by suDecripr R ^ the 

(6. 30) 

( 6=31 ) 

At mi dspan - denoted by eubecript ' m ' , allowance le made 
for a di Bcont 1 nui ty in elope of the circulation tunction- 
Note that the circulation t Line t ion is continuous ar 
midspan, but its derivative is di scont i nuous at this 



hiiJbPAN ANiJ wiNG j i h 



Hi one rne riqnT: wino 



D oun G a.r y c on g i r i on s are 



Tr^d = o 



AH 



= Pp':!) = o 



i ocat i on , 



F= BASIC Ar4D ADDITIOrlAL CIRCULATION 

An under st andi ng of the concepts of basi 
additional circulation is necessary in the analysis 
problem. These concepts are briefly reviewed here. 

Consider the circulation function [^(0,6) of a 
w?. no at some arbitrary angle of attack <=< , This 
expressed in terms of a “basic" circulation f 
and a "specific" additional circulation f 
as follows 

t ( cL - g) Pp, 



ano 



of thi s 



gi ven 
may be 



unct ] 



inci 



on 



( 6 



3 C 1 



- c< 






termed the additional angle of attack. 



where 



X =5 



a 3 HP 



It 


is useful to designate the 


product 


(C^ - ) |T ( 

t‘ H 




he 


'■‘net'* additional circulation 


and to 


d e 3 1 G rs a t e 


the 



G Lia n 1 1 1 V H- (0,6) it 33 i i 



t: P3 



aoGi -rj Qnai 



: 1 r c Lu a. n 1 



L-l 



3 rent i at i no 3qu.a.t i on ( 6 , 3‘u:i ) wi t. h P33p3c t to o< g a vb: 



r — L ^ = 6 ) J = O -f- ( 0 5 0 i 

doc 






Notice that [*^(0=6) hae the character of circulation 
per radian, and that it is nufrier i cal i y equal to the value 
which the net circulation would r^B.c:h if iineariy 
extrapolated to an additional angle 
radian. In fact, the additional angle 



oi* a'C'cacK of one 
of attack ( c< oC-v 



13 reeuricreG to email vaiuee; -Therefore^ tne net 
additional circulation remaine email in comparieon with the 
specific additional circulation. 

All other aerodynamic parametere which are proportional 
to circulation follow a linear law analogous to that stated 
in eq ua t i on e (6,32) an d ( 6 ,33), F or ex amp 1 e , the wing 1 i f r 
coefficient C* confer me to the exDreeeione 



C| = Cj_g 4- (o< - o<p)Ci £ 



B' '-LA 









= O -i~ L- s • 

Lf 



( 6 , 3o ) 
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note 



H Q o 1 n 

coet f i ci er 
(dC|_/di50 = 

e. doi ti one. 1 
addi t i onal 
.hoi-^eve.r" , 

ctec 
Basi c 



me speciric 



\na 



L'j - 15 idePiTiicai wion me lift curve sicoe 

value that the net 



J. ! ! 



u.'nericaiiy eouai to the 



lift would reach if 

angle of attack of one radian. Once again, 
tne aQQition ai angj 
? to 5 rn a 1 1 v' a i u e 5 • 
and additional 



ii. near IV ex or apoi acea oo ar 



attack “ ^w> 



44 



A . PRiJHLLri SPEU i F I CAT i ON 

As mentioned previously, the calculation sequence ma 



involve oeter mi ni nq the wing slope function 



II '—M ! ! 



specified basic lift di st r i but i on , or determining the I if 
distribution from a specified wing slope function. G 
interest nere are bo bn cne oasic and aooicionai iit 
d i s t r 1 bu 1 1 on s . 

The values on the right wing span are utilized in ai 
c a 1 c u 1 a X. 1 o n s • T n e r i g n t wing span is n i v* i d e g i n x o e i g p 
chordwise and eight spanwi se stations. 
siMxy— four control points over the semi— span. 



f = J. JL ! 



the condition of no flow t hr ouch 



wi n* 



!_= ! ! c*. •_ r=r s 



induced velocity must 
f r ee stream ve 1 oc i t y 



L_* o>. JL Cl ! j c: 



me normal comqonenx 



w. 



'cO 



= 



Setting the free stream velocit' 
approK i mat 1 ng sin<X with for sm< 

qi ves 



one un: 



ci ■_ c= u- r*. 






^ = W = W 

P P 



w„ 

5 -' 



(7.2 



By analogy with equation (6.32) 
function at each control ooint mav 



'eDresenxeD as 
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I! i 



IWQ COiT'OO 






w 1 t: n 



:ne oasic iii 



3.m 



"C r» a 3 Q Q 1 TI 1 Q r* 3 i i. 2 '1' TI ^ r“ 0 5 p a C T1 1 V a 1 y' 



S » I ! ^ 



^r-,. O \ Cl / 



(7,3 i 



where "i-iB the slope of 

accord! nq to the relation 



• inq cafTioer SLs.rtace 



Gi :4 



— 2 / K } tj 



wn lie 



.1 = — ’ X 

DH 



I / , Oi 



■efieccs Une L-niror 



3nae i r? Wj 



aoQiiiionai anpie 0 "i auTia(_5<ii 
Note that the right 



. ope 



o ^ o (_? a c 1 o r 






consists of the su.iTi of eoLiations (4.Z4) and (4,39), t 
solving for the additional lift^ the only unknowns 



eauati ons <4. 24' 



[4,39) ar 



I I 5 



e derivatives 



:ne 



c 1 r cui at 1 on t unction. I nese derivatives a 






re e-: pressed in 
the c 1 r c u. 1 a 1 1 on . 



!l! 



w 1 1 n SI X t y' to u. r Lt n r*i n g w n v' a i l? e s o i 
Q u a 1 1 o n s are required tor a s o 1 lj 1 1 o n 



Dener ates an 



c 1 r cLii at 1 on = s i x t y — -r our 

t-3.ch control poinit 
equation involving terms of 
d 1 s t r 1 b u 1 1 on . This p r od uc es 

equations in sixty— four unknowiis which musi 
SI mul taneousl y to determine the circulation Vc 
solution involves the inversion of a 64 bv 






ui at 1 on 



ne 


un known 


SI X 


ty-f ou?^ 


be 


sol vea 


ue s « 


This 


64 


m a. t r 1 X 3 



4 tj 



lii K Q r" ^ = 1 5 T or~ ti n ^ o (3. r~ i -3 1 ci ^ i“~ x v' 3. ii x v' 3 3 o t iz, n ^ cr x ci li x 3 n i c* n 

i* u. n cr ti 1 o r? x n n ^ rn = o "f ti n ^ p i k Ps o w n c i p q l? i s ^ x o n '•/ 3. x 33 -r^ r 0 



;e 0 oec 



are derivea eex 



i^jpen seiving the ha? 
D r e e e li re q x e t r x d u. t x on = 



J. * U «•! 



ps •- 



rns un Known 



DDxen^ witn a SDecitxed 



• !••_•_ JL 5 • 



! p X e sox u.Ti on 



'oxves fTU 



! ! i L.*. i V_ J 



me wino ex ope 



xcauxon ov a 



/I i-W V 



64 matrix^ but inversion ot the matrix 



equi r eo . 



UHOKhwIbb UERIVAixVbb OF C I KCULA I xUN 
In deveiopinq 



xnite Gxtterence eou 



.la t_ A L_« 



ne, a polynomial 



le paesBG tnrough a numoer of points in the particular area 
of interest. The poiynomiai used in each case must be of 



cient oroer 






me repuireo derivative? 



exist, and it must satisfy the boundary conditions. The 
finite difference equations utilized for the central 
difference formulas require a point on either side of the 
point 1 Pi question. rl e n c e = the c e n t r a x o x t f e r e n c e e o l* a t x o r? s 
are valid for all mesh points except those that lie 
adjacenc to the leading and trailing 



edqes. Fon 









+ Dce^'j 



Di f f erent i ati ng ^ 



^ ~ A^ + A^(^ + Get'll 
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EqLJatlons i / . h) and (/.7) satisfy the bcundary conditions 
for the additional lift. These are that P vanish and 



the leading edge where ^ = O. Applying equation 
(7«6) at two grid points adjacent to the leading edge= 
solving for Ai and A-.- and substituting into equation 
(7=7) , a difference expression for P^ in terms of P i and P 
is obtained as follows 



= <1 /A^ ) r -i- (1 /3 ) 



V / = H ; 



Therefore^ for any mesh point adjacent to 



leadino eoqe, 



the following subscript notation is used for 

( j , 1 ) = <1 / ) L r< j , i ) 4- (1 /3 ) P C j , 2 ) 1 



( 7 = 



A s 1 rn 1 1 a r ^;,n a 1 v' s i =? a o p 1 1 e □ o o c n e o a s i c 1 1 f ti r e s u i ti *— in c n e 
foil Qwi ng 

= (i/A^) c-<A^^/2)2r^^j_(j> -^6r<j,i)D < 7 = 10 ? 

The procedure for determining the difference formula 
applicable near the trailing edge requires defining a 
function for the circulation that is numerically equal to 
W at the trailing edge with first and second derivatives 
equal to zero= This requires beginning with a cubic 



n = r - 0)^ 



( 7 . i 1 ) 



Di f f erent 1 at i ng once, 
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= - ( 1 / 2 > A 3 (Tr - ) 2 



L'l f t ereni: 1 at 1 ng aga. 



^<i><k = ^A3<ir -<jS) 



\ / B 1 1 / 



iwo nr a. 1 1 I no 



solving tor » and sLibst i tut i ng into (7,12), 
eKprossion for Y in terms of r^(j,7)and P(j^S) 



n o ■? T 



a airte?^enc 



s oonaine 



oi 1 ow^ 



r^<j5S) = ( i /A0) L (— n-/ 1 3) n!j=7) + ( 3/ 1 3 ) P( j = B) 3 



V / » i ^ 



Hi so 



n. < j ) 



= (-1/26) r<j ,7) + (27/26) Pcj ,8) 



K / = 1 O 



For mesh points >’ through 7 in the chordwise dii^ection 
entral differences are used which results in the follow in 






(j,i) = 



(l/A(^) C (-1/2) r< j ,i-l) + (1/2)P( j ,i+l) 3 <7, 



The difference equations for the trailing edge and 1 
zsntra.l ooints are also valid for the basic lift or abler 



SPANwISE DERIVATIVES OF CIRCULATION 
Alono the right wino tio, the cir 



:ui am on 



however 5 the first derivative is finiteB These condition 
are precisely analogous to the boundary conditions a.lon 



the leading edge for the additional li 



difference ecuation developed above for 



Hence, the sam. 



n n ess o oun o ar 






ca =0 - »- or 



riant wine 



p0 ( 1 , 1 ) = ( 1 /AS ) r P( i . I ) -h (l/3)P(2=i)3 



.• 0 D t r a I 



or an c a +or?TiLiiaa ara ai=D tna aa.fna 



Tq ■' J ? 1 > 



= ( 1 /AS) L <-l/ 2 > r< j-l . 1 ) + ( 1 / 2 ) 1 ^- 



j- I “i / *7 i o % 



At mid span* 
iHa Si op a ot 



iiowance is maas tor a. ai sconti nui t v 



T u. n c t: I on 



man -Cl 



previously^ Again = a polynomial suitable to satisfy the 



bc^undary conditions is assume 



heretore* i e 1 



(r-Po) =H^(0-0p> + (i/2)Ap(0 - 0( 



< 7 , 1 9 '> 



and 



F 0 — H-| ■+■ Q 0c}.* 



< / = 



Satisfying equation (7»19) at j = 6 and j = /, and solving 



for Ai and Ao « elves the needed coefficients i\ 



3uat 1 ons 



( 7 = 1 V ) and ( 7 , 20 ) . Now , ^pp i y ^ ng equat i on < 7 - 1 • 

m 1 G s p a n ana equation v 7 ■ 2 at m 1 d s p a n and at j ~ g 1 v' 0 s 






^(i) = (3/S)r(6,i> - (5/4)r'(7,i) -i- (15/8)r(8,i) (7.21) 



= (i/Ae)nr(6,i) - 3r<7,i) + 2r(8,i)j 

Pq(s,i> = (1/A0) c (1/2) r<6,i )-2r<7,i >-^(3/2) n 



(8.1)3 (7.23) 



oo 



D, t-IELD ELEMENT AiND iMAGb. ELtnsi.NT 



Hsl.^I 1 


that 


the only d 


1 f f e 


rences in the va 


«- -• <-• r-<- -5- I— . 

• i •— J •- CP i • 


left 


s e rn 1 s- p a. n a 


re t 


he a 1 q eb r a i c s 1 q < 


and QT, Th 


e-.^^ef 0 


re. w_ ^ may 


be 


thoL=ght Of 3.5 


parts 3 one 


part 


con si st i n 


c? 0 


f the c 0 n t r i b u. 


p* a r c 1 c u 1 a. r 


f 1 el 


d point on 


the 


r 1 ght semi -span 


soec i f i ed 


con t.f^ 


0 1 point, a 


nd a 


not.her part cons 


contr 1 but 1 


on t 


hat the I 


mage 


I'leid point ’ 


semi —span 


makes 


on the 


same 


control point 


concept ap 


plies 


to Wp% t 


he 


tr ai 1 1 ng vor t 1 c 


alo.no the 


c r a. 1 1 


ing edge. 


Equa 


tions ( 4 = 24 ) and 


mod i f i ed s 


1 iqht 


Iv so that 


the 


effect of a fie 



e- 



on i:n0 i Bt t. 

a Tno 301710 



zn.B rioniL BBcni'-BpB.n a.nci r.no bt^bct. 01* izs 103.00 on tino ioft: 
3 00 i —span may be calcuiated si mui t aneoo.el y , using Bymm.Btry 
relations^ Let the SLibscript ’1 ' j^epresent the right 

semi —span and let subscript ' 2 ' represent the left 

3 e m i — sp a n - Ma k i ng the substitutions for the s y m m e t r i c 

elements and ccllecting terms , equations ( 4 , 24 ) and ( 4 , 39 ) 
rnay be rew.»^itten as 




5 5 

• I ' 1 




i_/ 4 ' 

4 TrlAR 




+ 



(n+np") 









•_j 1 



tl Q L». " T1 J. O r*! 5 K / a ^ ) s \ / • Jli 3. r* G ( / a il! .J a W h tf " 0 V -3, 1 LI 3. H 0 G 



0 a c n c on ti r Q 1 .o o i n t: o v ar n n 0 00 m i — Bp an = 
equationB in the siKtv— four unknown values 



r o G u. 00 =1 > *■ n V — -j- Q !J 



c 1 r o u. 1 a "c 1 on -r u.n c u 1 on , 



ine soiLimon 






equ.ations prodLices the desired val ues o-f the circLilation 
-fu.net ion at the mesh points over the right semi— span = 

Ea BASIC LIFT 

In the basic li-ft probi em a pressu.re d 1 st.^^ i bLit i on is 
specitied over the wing a leaving the corr espondi ng wing 
slope functio.n to be determined* In this case, eq’uation 
( 6 a 22 ), rearranged a gives the chordwise derivatives o^ the 
circulation -fLinction at each mesh point on the wing, as 
-f ol I ows 



-C (-^) Sm 



8 



\ / a ; 



At the leadi.ng edge* eqiuation (6-2*H) applies 



r0^L<-’> - f c(i) pjp 



t / a } 



Now, equations (7= 10) 5 (7*14)* and (7, 16) are e.-< pr essi ons 

for in terms c-f the Lin known circLiiation val Lies* These 

eoLiations are repeated belov^i with known valL»es col 1 ected on 
the left side of each eGLiation* 



= (i/A(z!i)6r( j , 1) 



( 7 a 1 O ) 



= = ( 1/A^) L (-1 /2) j , i“i ) -i- ( 1 /:.::) n< j « I -f- 1 ) J </,16: 

Pj^(j 58 ) = ( 1 /A^) C <-3/ 13) j , 7) (3/ 13) H( j , S) j (7.14^ 



H 1 i V' 3 i cr 5 on n r* ^ 



hand sides ot the above three 



i ) and ( / , 27 ) 



OT 



cn »-• y 



eoLiations are determinea ny equarions ». / 

The solution -for the circulation values requires inversion 
” ■ mat r i x , 

( / = 1 u ) = a n Q '. / » 1 o / are u. 1 1 1 1 z e g 
to calculate the spanwi se derivatives 
Ail variables 



' j o w ^ ^ 1 o r? 5 / a i / ) 



t unct 1 on : 



3P. the rioh-^ 



tne circLiiaziQ^ 



□ es G 



e Q u. a 1 1 o n s k / . 2 a ) and ( 7 a 2 t_j are now k n o w n » ! n e wine 



■ unct ion rna.'y' ce written as 



W_ \ K ^ ) 
D D 



- B(kp,k>r(k) 



i. / . 2 l 



Where y is a 64 by 64 matriK ot i 
whose e>^act detinition is implied by equations 
(7,25) a A more detailed anaivsis of matrix B 



1 uence coeri - 1 ci er 



<ef a 4 . 



Fa ADDITIONAL LIFT 

In the additional lift problem, 
pressure distribution over the wine 



= —1 a 



kjcr = J. 1 



oi 1 owi no o 



?r evi ousl y devel oped equat i ons 



are requirec 



r^(j,i) = < i /A^) c r< j , 1 ) + (i/3)r<j,2): 



<7. 9) 



n0(j,i) = <1/6^1) c <-i/2) To , 

,e> = ( i /L(^) c (-3/13) r< j 

Tq ( 1 5 i ) = (1 / A6 ) C 1 5 i ) -f- i. 
p0 ( j , i ) = (1 /A8 ) C < — 1 / 2 ) j — 
rg(S,i) = (1/A8) r (l/2)r(6,i 
r^ij) = (-1/26) Pi: j ,7) + <27 




! ne 



VV A • I y — - - - 



slope function in 



:his case mav be written 



I 



/ % r*f \ 

w-p,.-. I (r.. / 



/ » o4 ) 



w n e r e f-? 



w r I i_« = tz" e X a c 



iTi-arr 



OT 



Li ence 



1- 1 c i e 






on IS a Q a 1 n i rn d i i e b o v e q l? a s i o r 






and (7=25). Matrix A rnu.st be inverted in this case to 
deteroine the circulation function. 

Once the circulation function is known, its derivative 



:-pcr._ I 



1 n c o e n a •_ 1 o r 



' ^ 1 ■= f o u n c D V 
(7=9). (7,14). and (7.16), Eouation (6.26) 



IS then utilized to determine the pressure di str i but i on , 
Lj "C h e r a e r o G V n a m i c p a r a m e t e r s •may' t h e *n be r e a g i 1 'y c a i c 1 a c e g 
for u n e w i n g in *n l* e s 1 1 o n . 

■trices and calculazion procedures presence o in 



cne a.GQve oara^oraDns 



described in detail 



. n 



c.oncurrent analysis conducted by E. M. Barber LRef. 41 



VI 1 1 = SUhMARy 

Th^ •=-oiui:ion ot tn^ sv"=t.efr! or equa.iiions d^v^loped in 
this analysis is best SLtited +or processing by nacrix 
methods on a digital computer « A br ie-f summary o+ the 

p^■ocedural steps rBquirBd in the cal cul at i ons for the 

additional lift and the basic lift are given in tne 
^Qi lowing paragraphs, 

A= additional lift 

Biven Wp = -1; 

First, calculate matrix A utilizing equations ('7,9) , 
(7,14), (7-15), (7.16), (7-17), ( 7 , 1 S ) , an d (7,23) t og e t n er 

with the olantorm carameters in equations (7,24) and 



Next, invert matrix A and multiply bv vector to 

f t~ / t Q 

determine the circulation values- 

Nqw - determine the cnordwise derivative ot the 



circulation function utilizing equations (7.9), (7.14) 



, cri« = ^ 



(7r 16) , 

hinally, use equation (6,22) to calculate the pressure 
d i st r i bu t i on . 






Ba BASIC LIFT 

Gi V0.n P ( j = i ) • 

Firsts cs.lcuis.ts ths chord wi SB dsri v'stivs ot th^ 
circLilation -function using equation (6.22) , and the second 
d er i V a t i V e tact or u Bing eq u a t i on ( 6 = 2S ) = 

Next, compute the circulation function using equations 
(7=9) = (7=14)= and (7= 16) , 

Noi.-j , c omp u t e ma t r i x B , 

Finally, multiply matrix B by the circulation vector to 
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IX = CuNCLUSiUNH 



A1 thoLig 



-sr’e s.iTsO i & o 
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pnoaontod i 
tho QV0?‘-”a.i 1 



ad0' 



p e r f a r m a n c e 



h no numericai results are yet available^ 
rounds Tor con’ll dence in the validity ot 
method presented in this thesis, 
are based on sound mathe-mat 1 ca j procedures 
. This analysis^ along with the parallel e 
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has satisf 
effort * 
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